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Wood  combustion  experiments  were  carried  out  to  determine  the  effect  of  ignition  technique,  biomass  load  and 
cleavage,  as  well  as  secondary  air  supply,  on  carbon  monoxide  (CO),  total  hydrocarbon  (THC),  particulate  matter 
(PM10)  and  particle  number  emissions  from  a  woodstove.  Wood  from  two  typical  tree  species  in  the  Iberian  Pen¬ 
insula  was  selected:  pine  ( Pinus  pinaster )  and  beech  ( Fagus  sylvatica).  The  highest  CO  and  total  hydrocarbon 
emission  factors  (EFs)  were  observed,  respectively,  for  pine  and  beech,  for  high  and  low  fuel  loads.  The  highest 
PM10  EF  was  recorded  for  the  operation  with  low  loads  for  both  woods.  Secondary  air  supply  produced  the  lowest 
PM10  emission  factors.  The  top  ignition  can  decrease  the  PMi0  EF  to  less  than  half  when  compared  with  the  com¬ 
mon  technique  of  lighting  from  the  bottom.  The  lowest  particle  number  emission  factors  were  observed  when 
operating  with  high  loads  of  split  beech  logs  and  when  using  secondary  air  supply  during  the  combustion  of 
pine.  Regarding  particle  number  distributions,  the  highest  geometric  mean  diameter  (Dg),  for  both  woods, 
were  observed  when  operating  with  high  loads  (with  split  and  non-split  wood). 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Currently,  most  developed  countries  are  heavily  dependent  on  fossil 
fuels.  In  this  context,  policies  that  encourage  the  use  of  biomass  as  a  re¬ 
newable  energy  source  are  promoted  [1,2],  Despite  the  recognised 
advantages  of  the  use  of  biomass,  residential  wood  combustion  for 
heat  production  has  been  assessed  to  be  a  major  source  of  fine  particle 
mass  emissions  [3,4]  especially  during  wintertime  [5,6],  The  problems 
associated  with  this  sector  occur  particularly  during  winter,  not  only 
for  being  the  period  in  which  residential  wood  combustion  is  most 
commonly  used,  but  also  due  to  weather  conditions  of  stagnation. 
Thus,  residential  wood  combustion  can  result  in  local  particulate  matter 
levels  comparable  to  those  registered  in  high  traffic  areas  [7],  Emissions 
from  residential  wood  combustion  are  generated  only  a  few  metres 
above  the  ground;  consequently,  these  emissions  have  not  enough 
time  to  dilute,  chemically  oxidise  and  react  before  the  population  is 
exposed  to  them  [8], 

In  Portugal,  it  is  estimated  that  2  Mton  of  wood  is  annually  used  in 
residential  combustion,  which  contributes  to  particulate  emissions 
representing  around  30%  of  the  total  PM10  emitted  by  the  diverse  sec¬ 
tors  of  activity  in  the  country  [4j. 

Airborne  particles  arising  from  residential  combustion  have  harmful 
effects  on  public  health  [9-11].  Some  studies  have  been  conducted  to 
evaluate  the  toxicity  and  mutagenicity  of  particulate  matter  resulting 
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from  domestic  biomass  combustion  [12-14].  Epidemiologic  studies 
have  linked  fine  particulate  air  pollution  with  negative  health  effects 
like  respiratory  and  heart  diseases  or  even  premature  dead  [15,16].  Sev¬ 
eral  studies  have  shown  that  the  submicrometer  sized  particles  (<  1  pm) 
are  predominant  in  PM  emissions  from  this  source  [17],  The  smaller  the 
size  of  the  particle  the  greater  can  be  the  risk  for  human  health.  Huang 
et  al.  [18]  have  been  following  cardiovascular  disease  patients  during 
two  years  and  observed  a  reduction  in  heart  rate  variability  index  due 
to  reduced  exposure  to  PM2.5  and  black  carbon. 

Taking  into  account  the  impact  of  these  emissions  and  the  need  for 
compliance  with  legal  norms,  a  rigorous  quantification  and  characteri¬ 
sation  of  emissions  from  this  sector  are  necessary.  Additionally,  and 
since  previous  works  have  suggested  that  wood  and  the  type  of  com¬ 
bustion  appliance  exert  a  major  influence  on  emissions  from  this  sector 
[19-24],  it  is  necessary  to  obtain  values  for  Southern  Europe  in  order  to 
use  more  representative  and  specific  locally  acquired  data  on  emission 
inventories  and  source  apportionment,  photochemistry  and  climate 
change  models. 

The  technology  used  in  domestic  biomass  combustion  has  a  prepon¬ 
derant  influence  on  the  observed  emission  profiles.  This  is  especially  evi¬ 
dent  when  emissions  from  manually  and  automatically  operated  systems 
are  compared.  Schmidl  et  al.  [24]  reported  high  variability  between  man¬ 
ual  and  automatic  systems.  PM10  concentrations  obtained  by  these  au¬ 
thors  were  between  12  and  21  mg  m“3  for  automatic  systems  (pellet 
boilers  and  pellet  stoves )  and  in  the  range  1 1 1  -1 5 1  mg  m~  3  for  manually 
operated  systems  (typical  stove  and  stove  with  secondary  air  intake). 
Johansson  et  al.  [25]  reported  emissions  50  times  higher  in  stoves  and 
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fireplaces  compared  to  more  controlled  equipment.  Johansson  et  al.  [23] 
observed  concentrations  of  PM  in  old  boilers  about  1 80  times  higher 
than  those  emitted  by  modern  pellet  burners.  Schmidl  et  al.  [24]  also 
found  variability  in  the  composition  of  the  emitted  particles.  For  manually 
operated  equipment  the  authors  showed  that  70%  of  the  PM10  mass 
consisted  of  total  carbon  (TC).  For  the  automatic  combustion  equipment, 
the  emitted  particles  consisted  mostly  of  potassium  inorganic  salts.  Fine 
et  al.  [26]  reported  that  40  to  100%  of  the  fine  particle  mass  was  organic 
material  when  a  conventional  fireplace  was  used. 

Although  the  emissions  from  residential  wood  combustion  have 
been  reported  in  previous  studies,  there  is  still  a  limited  understanding 
of  the  influence  of  operating  practices  on  these  emissions.  They  vary 
considerably  among  studies,  affecting  the  amount  and  composition  of 
the  flue  gas  emissions.  In  addition  to  batch  operation  mode  and  log 
size,  as  well  as  fuel  quality  and  air  supply  strategies,  it  is  also  necessary 
to  take  into  account  that,  when  studying  the  effect  of  practices  on  emis¬ 
sions,  each  appliance  has  its  optimal  best  operating  conditions.  Different 
operating  practices  are  also  relevant  with  respect  to  the  effects  on 
human  health.  Leskinen  et  al.  [27]  studied  the  toxicological  responses 
of  the  fine  particles  resulting  from  three  different  combustion  condi¬ 
tions,  namely  efficient  combustion,  intermediate  combustion,  which 
simulated  improper  air  staging  and  partial  load,  and  smouldering  com¬ 
bustion.  The  different  combustion  efficiencies  were  simulated  in  a  wood 
chip  grate  combustion  furnace.  In  efficient  combustion  conditions,  the 
emitted  fine  particles  had  the  highest  potential  to  cause  cell  death. 

The  objectives  of  the  work  presented  here  were  to  identify  the  influ¬ 
ence  of  operation  and  real-life  conditions  on  emissions  from  residential 
wood  combustion.  The  investigated  parameters  were:  (i)  ignition  tech¬ 
nique  (upside-down  and  bottom-up  lighting),  (ii)  hot  start  versus  cold 
start,  (iii)  fuel  load  and,  in  the  case  of  high  load,  degree  of  cleavage  of 
the  logs,  and  (iv)  secondary  combustion  air  supply.  This  evaluation  al¬ 
lows  assessing  which  practices  offer  the  greatest  emission  reductions 
and,  in  consequence,  could  constitute  an  important  step  towards  an  im¬ 
portant  reduction  of  these  emissions  in  South  Europe.  Moreover,  it  con¬ 
stitutes  an  opportunity  to  characterise  more  thoroughly  the  emissions 
from  residential  wood  combustion  by  replicating  operating  conditions 
commonly  carried  out  in  southern  European  countries  and  to  point 
out  measures  that  can  reduce  the  emissions  from  this  sector. 

2.  Methods 

2.1.  Combustion  facility 

The  combustion  tests  were  performed  with  a  typical  Portuguese 
stove  (power  output  9.6-18.2  kW),  commonly  used  for  domestic 
heating,  manufactured  by  Solzaima  (Sahara  model).  The  rated  output 
of  the  stove  is  14.0  kW.  This  equipment  is  made  of  stainless  steel  with 
a  cast  iron  grate  (Fig.  1).  The  combustion  chamber  has  a  volume  of 
0.093  m3,  corresponding  to  a  height  of  0.44  m,  a  width  of  0.59  m  and 
a  length  of  0.36  m.  It  is  a  batch  operated  furnace  with  manual  control 
of  combustion  air;  besides  performing  heating  by  radiation  and  natural 
convection,  it  has  also  been  designed  to  include  air  heating  by  forced 
convection. 

The  combustion  air  enters  the  combustion  chamber  through  a  regu¬ 
lating  device  located  at  the  bottom  of  the  equipment,  in  the  ash  collec¬ 
tion  hopper,  below  the  grate  of  the  stove.  The  flow  rate  of  primary  air 
entering  in  the  combustion  chamber  was  monitored  by  a  mass  flow 
meter  (Kurz  Model:  500-40  0.0  P-2).  This  model  was  designed  with 
built-in  inlet  piping  with  the  same  cross  section  of  the  stove's  air  intake 
to  ensure  underneath  air  flow  identical  to  that  originally  fed  in  under 
the  grate.  The  driving  force  for  the  air  flow  rate  throughout  the  combus¬ 
tion  chamber  is  the  natural  convection  resulting  from  the  up-flowing 
stream  of  combustion  hot  flue  gases  throughout  the  chimney.  The  nat¬ 
ural  draught  mainly  depends  on  flue  gas  temperature,  mass  flow  and 
ambient  temperature.  This  fact  may  cause  some  variability  on  the  re¬ 
sults.  It  should  be  noted  that  the  woodstove  was  operated  under  similar 


conditions  to  those  practiced  by  household  users.  Contrarily  to  what 
happens  in  Portugal,  some  other  countries  adopted  regulations,  which 
indicate  that  the  chimney  must  be  able  to  provide  a  draught  of  at  least 
1.2  mm  water  gauge  (12  Pa).  The  stove  is  equipped  with  a  vertical  chim¬ 
ney  with  an  internal  diameter  of  0.2  m  and  4.0  m  tall,  insulated  with 
ceramic  fibber  (Cerablanket,  128  kg  m-3,  0.025  m,  0.23  W  m1  K-1)  in 
the  first  1  m  above  the  exit  of  the  combustion  chamber  in  order  to  guar¬ 
antee  operational  safety. 

The  stove  was  loaded  manually  by  putting  each  wood  batch  over  the 
fixed  grate  installed  in  the  combustion  chamber.  The  grate  was  placed 
over  a  weight  sensor  (DSEUROPE  Model  535QD-A5),  which  allowed 
continuous  monitoring  of  the  fuel  mass  during  the  combustion  runs. 

The  combustion  flue  gas  temperature  was  monitored  using  K-type 
thermocouples  at  several  locations  along  the  system  (combustion 
chamber  and  chimney).  A  more  detailed  description  of  the  combustion 
facility  and  appliance  can  be  found  elsewhere  [28]. 

2.2.  Combustion  procedures  and  conditions 

The  combustion  experiments  were  planned  in  order  to  evaluate  the 
influence  of  the  following  operating  variables:  ignition  technique  (up¬ 
side-down  and  bottom-up  lighting),  fuel  load,  in  the  case  of  high  load, 
split  (S)  and  non-split  (NS)  logs,  and  secondary  combustion  air  supply. 
To  investigate  the  influence  of  fuel  load  on  flue  gas  emissions  three 
different  conditions  were  tested:  low,  medium  and  high  load,  according 
to  the  fuel  weight  (Table  1).  These  experiments  were  initiated  by  put¬ 
ting  a  batch  of  fuel  over  the  bed  of  hot  char  from  a  previous  combustion 
experiment,  when  the  temperature  in  the  combustion  chamber  was 
around  100  °C  (±20  °C).  All  the  tests  to  evaluate  the  effect  of  load  and 
log  cleavage  (S  and  NS  logs)  on  emissions  correspond  to  a  condition  of 
“hot  start”.  The  wood  was  cut  into  logs  of  30-40  cm  in  length,  whether 
for  split  or  non-split  logs.  Taking  into  account  that  the  log  size  influences 
the  combustion  quality,  the  wood  was  purchased  from  a  local  supplier 
with  a  large  market  share  in  an  attempt  of  replicating  the  Portuguese 
householder’s  practices.  The  duration  of  each  combustion  cycle  depended 
on  the  operating  condition,  load  and  wood  type,  ranging  from  30  to 
95  min.  The  fuel  ignition  techniques  included  the  conventional  ignition 
from  the  bottom  and  ignition  from  the  top.  To  initiate  the  combustion 
experiments  using  the  bottom-up  method  of  lighting,  two  pine-cones 
were  kindled  with  a  lighter  and  the  fuel  load  was  put  on  the  veiy  top. 
The  upside-down  fire-starting  method  was  achieved  using  small  pieces 
cut  from  the  same  wood  being  burned  and  pine-cones  cracked  on  the 
top  of  the  batch  of  logs.  These  experiments  were  made  with  wood  batches 
of  2  kg  for  both  fuels  (medium  load).  The  total  number  of  batches  for  test¬ 
ing  each  ignition  technique  was  between  three  and  five  (Table  1 ). 

To  investigate  the  influence  of  secondary  combustion  air  supply, 
experiments  were  performed  using  a  secondary  air  inlet  (3.8  cm  inter¬ 
nal  diameter)  at  the  axial  axis  of  the  rear  wall  of  the  combustion  cham¬ 
ber  and  located  at  16  cm  above  the  grate  of  the  furnace,  in  order  to 
promote  a  secondary  combustion  zone.  The  position  and  dimension  of 
the  secondary  air  supply  was  defined  after  discussion  with  the  manufac¬ 
turer.  Secondary  combustion  air  entered  the  combustion  chamber  by 
natural  convection  as  a  result  of  the  chimney  drag  force  induced  by 
the  up-flowing  hot  combustion  flue  gases.  The  aim  was  to  enhance 
the  mixing  between  pyrolysis  gases  from  the  thermal  decomposition 
of  the  wood  and  combustion  air,  in  order  to  achieve  improved  burnout. 

2.3.  Wood  fuel  characteristics 

Biofuels  were  selected  from  tree  species  typical  of  the  Iberian  Penin¬ 
sula  and  commonly  used  in  households  for  domestic  heating:  Pinus 
pinaster  (softwood)  and  Fagus  sylvatica  (hardwood).  The  moisture  and 
ash  content,  together  with  the  elemental  composition,  are  shown  in 
Table  2.  The  determination  of  ash  and  moisture  contents  was  performed 
in  accordance  with  CEN/TS  14775  and  CEN  14774-3,  respectively. 
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Fig.  1.  Schematic  representation  of  the  experimental  installation.  Dashed  line  —  Electric  circuit,  continuous  line  —  Pneumatic  circuit.  A  —  Stove;  B  —  Crate  of  the  stove;  C  —  Load  cell  (weight 
sensor);  D  —  Air  flow  meter;  E  —  Thermal  insulation  of  the  exhaust  duct;  F  —  Exhaust  duct  (Chimney);  G  —  Water-cooled  gas  sampling  probe;  H,IJ,K,L,M,N,0,P,Q,R.S  —  Gas  sampling  and 
analysis  system,  and  control  and  data  acquisition  system;  T  —  Dilution  tunnel;  U  —  Aerosol  size  distribution  probe  (PCASP-X);  V  —  Venturi  system;  W  —  Pitot  tube;  X,Y  —  TECORA  PM , f) 
sampling  system;  Z  —  Fan. 


2.4.  Flue  gas  monitoring 

The  composition  of  the  flue  gas  was  monitored  using  a  system  that  in¬ 
cludes  a  water-cooled  sampling  probe,  a  set  of  gas  conditioning  and  dis¬ 
tribution  units,  and  on-line  gas  analysers.  A  sample  of  the  flue  gas  was 
dried  by  condensing  the  water  vapour  and  other  gaseous  condensable 


compounds  in  an  impinger  immersed  in  an  ice  bath.  Then,  the  dry  gas 
was  distributed  to  a  set  of  on-line  gas  analysers.  The  continuous 
determination  of  the  02  concentrations  was  made  by  a  paramagnetic 
analyser  (ADC,  02-700  model,  with  a  Servomex  module).  C02  and  CO 
were  monitored  by  a  non-dispersive  infrared  analyser  (Environnement 
MIR  9000  model).  For  monitoring  total  hydrocarbons  (THC,  expressed 
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Table  1 

Details  of  the  combustion  experiments. 


Cleavage 

Load 

Fuel  burnt 

Number  of 

Flue  gas  temperature  (°C)  in  the 

Flue  gas  temperature  (°C)  at  the 

Rate  of  fuel  conversion 

(kg) 

batches 

combustion  chamber 

chimney  exit 

(kgh-1) 

Pinus  pinaster 

Split 

Low  load 

1 .0-1.2 

5 

284  ±  41 

135  ±  20 

1.22  ±  0.17 

Medium  load 

1.9-2.1 

5 

358  ±  36 

164  ±  12 

1.87  ±  0.38 

High  load 

3.7-4. 1 

3 

417  ±  36 

205  ±  18 

3.30  ±  0.55 

Secondary  combustion  air 

2.2-2.8 

5 

421  ±  51 

366  ±  26 

2.51  ±  0.42 

Bottom-up  ignition 

1. 9-2.1 

4 

329  ±  48 

151  ±  18 

2.25  ±  0.18 

Top-down  ignition 

2.1-2.2 

3 

356  ±  18 

163  ±  7 

2.31  ±  0.20 

Non-split 

High  load 

3.9-4.0 

3 

393  ±  83 

180  ±  18 

2.43  ±  0.05 

Fagus  sylvatica 

Split 

Low  load 

1 .0-1.2 

3 

252  ±  9 

113  ±  10 

1.12  ±  0.06 

Medium  load 

1.8-2.1 

3 

412  ±  27 

185  ±  2 

2.38  ±  0.07 

High  load 

3.8-3.9 

3 

440  ±  51 

216  ±  13 

3.32  ±  0.52 

Bottom-up  ignition 

1. 8-2.0 

4 

362  ±  42 

167  ±  17 

2.78  ±  0.13 

Top-down  ignition 

1 .8—1.9 

3 

350  ±  36 

149  ±  10 

2.13  ±  0.34 

Non-split 

High  load 

3.8-3.9 

3 

413  ±  2 

202  ±  8 

2.80  ±  0.09 

as  CH4),  a  heated  probe  and  sampling  line  (at  180  °C)  were  used  together 
with  an  analyser  with  flame  ionisation  detector  (Dyna-FID,  model  SE- 
310). 

2.5.  Particle  measurement 

Sampling  of  particulate  matter  with  aerodynamic  diameters  below 
10  pm  (PM10)  was  performed  in  a  dilution  tunnel  installed  downstream 
of  the  chimney  under  isokinetic  conditions.  The  entrance  of  the  flue  gas 
into  the  dilution  tunnel  is  made  by  suction  and,  for  this  purpose,  the 
input  of  the  dilution  tunnel  is  wider  than  the  output  of  the  chimney. 
This  configuration  simultaneously  promotes  dilution  and  effective 
mixing  due  to  the  turbulence  generated  [24].  The  dilution  of  the  flue 
gas  to  sample  the  exhaust  particles  is  widely  used  [17,21,29-33]  since 
it  allows  the  simulation  of  the  rapid  cooling  and  mixing  that  occurs 
when  the  exhaust  gases  are  released  into  the  atmosphere.  With  this 
aim,  the  inlet  of  the  dilution  air  (ambient  air)  takes  place  at  the  exit  of 
the  chimney.  The  dilution  tunnel  has  0.20  m  internal  diameter,  and  a 
length  of  1 1  m.  The  particulate  matter  sampling  point  was  located  at 
10  m  downstream  the  dilution  tunnel  entrance.  The  sampling  train  in¬ 
cluded  a  specific  sampling  head  (PM10),  a  pump,  and  a  control  and 
data  acquisition  system,  all  part  of  a  TCRTECORA  (model  2.004.01)  in¬ 
strument.  The  equipment  has  been  operated  at  a  flow  of  2.3  m3  h-1 
onto  quartz  filters  (47  mm  diameter).  To  ensure  sampling  under 
isokinetic  conditions  a  stainless  steel  probe  with  an  internal  diameter 
of  1 1  mm  was  used.  Two  consecutive  filters  were  collected  during 
each  combustion  test.  Since  the  replacement  took  less  than  1  min,  the 
overall  particle  emission  factor  for  each  combustion  cycle  was  consid¬ 
ered  not  to  be  significantly  affected.  The  temperature  in  the  particle 
sampling  point  was  in  the  range  25-35  °C,  and  the  relative  pressure 
was  typically  at  98.1  Pa  below  the  atmospheric  pressure.  The  volumetric 
gas  flow  rate  throughout  the  tunnel,  and  respective  combustion  gas  di¬ 
lution  ratio  were  calculated  from  the  mean  gas  velocity  in  the  cross  sec¬ 
tion  of  the  dilution  tunnel.  This  mean  gas  velocity  was  estimated  from 


Table  2 

Moisture,  ash  and  elemental  composition  of  the  wood  used  as  fuel  in  the  combustion 
experiments. 


Pinus  pinaster 

Fagus  sylvatica 

(wt.%,  as  received) 

Moisture 

9.90 

9.60 

(wt.%,  dry  basis) 

Ash 

0.40 

1.80 

C 

51.40 

47.97 

H 

6.20 

6.26 

N 

0.16 

0.04 

S 

bdl* 

bdl* 

O" 

41.84 

43.93 

*below  detection  limit  of  0.01  wt.%;  **by  difference. 


the  differential  pressure  monitored  by  a  Pitot  tube  and  respective  pres¬ 
sure  sensor  (Testo  AG  808)  and  a  K-type  thermocouple  [29,34].  The  gas 
flow  velocity  was  6.3  ±0.1  ms-1  and  the  dilution  ratios  applied  to  the 
flue  gases  from  the  woodstove  were  around  25:1.  The  volumetric  flow 
ratio  was  obtained  by  multiplying  the  flue  gas  velocity  by  the  cross  sec¬ 
tion  of  the  dilution  tunnel.  The  dilution  factor  was  estimated  by  dividing 
the  flow  rate  in  the  dilution  tunnel  by  the  flow  of  exhaust  gases. 

After  each  combustion  experiment,  the  PMi0  sampling  head  was 
cleaned  to  minimise  contamination  of  the  subsequent  experiments. 
Background  samples  were  collected  to  assess  any  contamination  arising 
from  the  dilution  air.  The  average  of  several  blanks  was  subtracted  from 
measured  values. 

Gravimetric  analysis  was  performed  with  a  microbalance  (RADWAG 
5/2Y/F  with  an  accuracy  of  1  pg  after  48  h  equilibration  in  a  room  with 
controlled  temperature  and  humidity.  Quartz  fibre  filters  used  in  this 
work  were  ignited  at  500  °C  for  6  h  to  eliminate  organic  contaminants 
before  sampling.  The  baked  filters  require  a  more  careful  storage,  as  eas¬ 
ily  absorb  gaseous  organic  compounds  and/or  carbon  dioxide  from  the 
environment  in  substantial  quantities  [35],  Thus,  in  order  to  avoid  con¬ 
tact  with  sources  of  contamination,  exposure  to  ambient  air  after  calci¬ 
nation  was  avoided,  keeping  the  filters  in  petri  dishes.  The  transport  and 
storage  were  carried  out  in  a  desiccator.  Filter  weight  was  obtained  from 
the  average  of  six  measurements,  when  the  variations  were  less  than 
0.02%. 

The  particle  number  concentrations  and  particle  number  size  distri¬ 
butions  are  also  important  parameters  to  develop  a  quantitative  assess¬ 
ment  of  their  impacts  on  both  human  health  and  the  environment. 
The  particle  number  size  distributions  were  measured  with  a  laser 
spectrometer  (Passive  Cavity  Aerosol  Spectrometer  Probe,  PMS  Model 
PCASP-X)  in  a  second  dilution  tunnel  with  an  internal  diameter  of 
0.07  m  and  1.56  m  length.  The  device  measures  the  size  distribution  of 
latex  particles  with  nominal  optical  diameters  between  0.1  and  10  pm 
in  31  discrete  channels.  To  obtain  the  exact  number  concentration  of 
particles  in  each  channel  it  is  necessary  to  made  several  corrections  to 
the  spectrometer  counts  [36].  Since  the  optical  particle  counter  was  cal¬ 
ibrated  with  latex  particles  with  a  refractive  index  of  1.58-0.0i,  it  is  nec¬ 
essary  to  make  the  refractive  index  correction  for  the  aerosol  sampled. 
The  correction  of  the  raw  size  bins  in  order  to  obtain  the  exact  number 
concentration  versus  exact  size  was  done  using  a  programme  based  on 
the  Mie  Theory.  The  refractive  index  is  a  parameter  used  to  describe  the 
scattering  and  absorption  characteristics  of  a  particle.  Different  method¬ 
ologies  have  been  proposed  to  estimate  the  particle  refractive  index 
[37-39],  Due  to  the  lack  of  information  on  optical  properties  of  the  sam¬ 
pled  particles,  the  methodology  proposed  by  Levin  et  al.  [39]  seemed  to 
be  the  most  suitable,  since  it  was  intended  to  estimate  the  refractive 
index  of  particles  from  their  chemical  composition,  assuming  that  the 
PMio  constituents  are  present  in  the  form  of  certain  chemical  com¬ 
pounds,  having  a  specific  density  and  typical  refractive  index.  Since  it 
was  found  that  carbonaceous  constituents  accounted  for  the  largest 
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particle  mass  fraction  (42-96%),  the  refractive  index  was  estimated 
based  on  these  compounds.  It  was  also  necessary  to  make  the  adjust¬ 
ment  of  the  sample  volume  according  to  the  altitude  of  the  sampling 
point.  This  correction  was  made  as  recommended  by  the  manufacturer. 

3.  Results  and  discussion 

3.1.  Combustion  conditions 

For  each  combustion  test,  the  continuous  monitoring  of  the  follow¬ 
ing  parameters  was  performed:  1 )  fuel  load  on  the  grate,  2)  combustion 
air  flow  rate  entering  the  combustion  chamber,  3)  temperature  of  the 
exhaust  gas  in  the  combustion  chamber  and  in  the  chimney,  and 
4)  02,  C02,  CO  and  THC  concentrations. 

The  combustion  cycles  were  characterised  by  three  main  stages, 
namely  ignition,  flaming  and  smouldering.  The  first  stage  corresponds 
to  the  heating  of  the  biomass  followed  by  drying  and  initial  steps  of 
devolatilisation,  without  the  existence  of  visible  flame.  This  initial 
stage  of  combustion  is  a  critical  step,  since  the  fast  devolatilisation  of 
the  fuel  originates  high  amounts  of  CO  and  THC  in  a  short  time.  These 
compounds  are  not  oxidised,  despite  the  high  concentration  of  02  avail¬ 
able.  This  can  be  explained  by  the  low  temperatures  in  the  combustion 
chamber  and  poor  mixing  between  the  combustion  air  and  flue  gas.  As 
long  as  the  02  concentrations  are  too  high  and  the  temperatures  in  the 
combustion  chamber  are  too  low  to  achieve  appropriate  burnout  condi¬ 
tions,  high  CO  and  THC  concentrations  in  the  flue  gas  are  detected.  The 
flaming  stage  is  characterised  by  de  volatilisation  and  combustion  of  vol¬ 
atiles  and  char  during  which  a  vigorous  flame  is  observed.  In  this  stage 
the  02  availability  in  the  flue  gas  and  sufficiently  high  temperatures  in 
the  combustion  chamber  provide  acceptable  burnout  conditions.  Dur¬ 
ing  this  phase,  the  oxygen  concentration  in  the  flue  gases  decreases  rap¬ 
idly,  as  well  as  the  THC  and  CO  concentrations.  On  the  other  hand,  the 
C02  concentration  in  the  flue  gas  increases  due  to  the  oxidation  of  vola¬ 
tiles  (Fig.  2).  At  the  end  of  the  combustion  cycle,  char  combustion  takes 
place.  In  the  smouldering  stage,  the  02  concentrations  in  the  flue  gas 
start  to  increase  and  the  temperatures  in  the  combustion  chamber 
decrease  again.  Consequently,  the  burnout  quality  decreases  and 
the  CO  emissions  increase.  Small-scale  combustion  appliances  use  natu¬ 
ral  draught.  Therefore,  in  the  ignition  phase  and  initial  steps  of 
devolatilisation,  large  amounts  of  CO  and  HC  are  formed,  but  the 
draught  of  air  is  scarce.  During  the  main  combustion  stage,  when  the 
draught  of  air  is  large  enough,  the  devolatilised  gases  are  oxidised.  If 
the  draught  is  too  strong,  i.e.  if  the  combustion  rate  is  too  high,  the 
devolatilised  gases  could  go  directly  through  the  chimney  without 
reacting,  due  to  insufficient  residence  time  in  the  combustion  chamber 
[40  j.  Thus,  two  periods  of  high  CO  and  THC  levels  were  identified  in  the 


present  work  and  in  other  studies  [28]:  a  first  peak  related  to  the  fuel 
fast  devolatilisation  after  reloading  the  combustion  chamber  and  a  sec¬ 
ond  peak  related  to  the  char  combustion  at  relatively  low  temperatures 
(<500  °C). 

The  maximum  CO  concentrations  for  the  combustion  of  pine  occurred 
during  the  operation  with  high  fuel  load  using  split  (7546  ppmv,  dry 
gases)  and  non-split  logs  (7661  ppmv,  dry  gases).  The  combustion  of 
beech  showed  the  same  trend,  also  with  high  values  during  the  operation 
with  high  loads  of  both  split  (8045  ppmv,  dry  gases)  and  non-split  logs 
(741 1  ppmv,  dry  gases).  The  maximum  THC  concentrations  were  also  re¬ 
corded  under  the  same  operating  conditions:  3401  ppmv  (wet  gases)  for 
pine  and  4031  ppmv  (wet  gases)  for  beech.  The  fuel  load  had  a  clear  effect 
on  emissions.  The  use  of  high  loads,  with  and  without  cleavage,  may  lead 
to  incomplete  combustion  due  to  the  insufficient  combustion  air  supply 
and  poor  mixing  between  the  combustion  air  and  flue  gas. 

Regarding  the  evolution  of  the  combustion  air  flow  admitted  to  the 
combustion  chamber  during  a  combustion  cycle,  there  is  a  sudden  in¬ 
crease  during  the  first  few  minutes  (5-10  min),  until  it  reaches  its  max¬ 
imum  value  and  then  decreases  gradually  until  the  end  of  the  test 
(Figs.  3  and  4). 

When  air  staging  was  tested,  it  was  observed  that  the  primary  com¬ 
bustion  air  flow  was  the  lowest  among  all  the  tests  performed.  The  flow 
rate  maximum  value  ranged  between  32  and  34  Nm3  h_1.  The  air  staging 
enables  to  reduce  the  combustion  air  supplied  directly  to  the  fuel  bed, 
i.e.,  the  flow  of  primary  combustion  air.  Reducing  the  primary  combustion 
air  allows  to  keep  the  temperatures  in  the  combustion  chamber  on  a 
moderate  level.  Maintaining  the  same  overall  combustion  air  ratio  allows 
to  use  reduced  amounts  of  primary  combustion  air,  which  results  in 
higher  amounts  of  secondary  combustion  air  entering  the  combustion 
chamber  and,  therefore,  higher  turbulence  and  an  improved  mixing  [41]. 

With  respect  to  biomass  consumption,  it  was  possible  to  distin¬ 
guish  two  phases.  In  the  first  one,  a  high  fuel  consumption  rate  was 
observed,  which  graphically  is  consonant  with  an  accentuated 
slope  (Figs.  3  and  4).  In  the  second  phase  a  slow  loss  of  mass  over 
time  was  registered,  which  corresponds  to  the  smouldering  stage 
of  the  combustion  cycle.  For  medium  load  combustion  conditions, 
it  seems  that  the  first  phase  occurred  during  the  first  30  min  and 
the  second  phase  from  this  moment  until  the  end  of  the  combustion 
cycle.  The  duration  of  both  phases  of  fuel  consumption  is  related  to 
the  duration  of  the  entire  combustion  cycle,  which  in  turn  is  a  func¬ 
tion  of  the  fuel  loading.  A  direct  relationship  between  temperature  in 
the  combustion  chamber  and  combustion  air  flow  rate  was  observed. 
When  the  air  flow  rate  entering  the  combustion  chamber  was  higher, 
the  duration  of  the  combustion  cycle  was  shorter.  The  highest  tem¬ 
peratures  were  reached  under  the  operation  with  high  loads,  whilst 
the  lowest  temperatures  were  observed  for  the  operation  with  low 
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Fig.  2.  Evolution  of  the  02,  C02  and  CO  concentration  in  the  flue  gas  at  the  chimney  exit 
during  the  combustion  of  a  batch  of  fuel  operating  with  medium  load  and  using  wood 
from  pine  as  fuel. 


Fig.  3.  Evolution  of  the  fuel  weight  in  the  grate  and  the  air  flow  rate  entering  the  combus¬ 
tion  chamber  during  the  combustion  of  a  batch  of  fuel  operating  with  high  load  (S)  using 
wood  from  pine  as  fuel. 
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Fig.  4.  Evolution  of  the  fuel  weight  in  the  grate  and  the  air  flow  rate  entering  the  combustion 
chamber  during  the  combustion  of  a  batch  of  fuel  operating  with  low  load  using  pine  as  fuel. 

loads  (Table  1 ).  The  highest  fuel  rate  conversion  was  observed  for 
the  operation  with  high  load  (S)  for  both  fuels,  whilst  the  lowest 
was  recorded  for  operation  with  low  load  (Table  1).  For  the  opera¬ 
tion  with  medium  loads  (with  and  without  secondary  air  supply),  the 
combustion  cycles  lasted  approximately  60  min.  Contrary  to  the  condi¬ 
tions  of  medium  loads,  for  which  the  duration  of  a  cycle  is  very  similar 
between  batches,  the  combustion  cycles  with  low  and  high  fuel  loads 
showed  highly  variable  durations.  This  variability  results  in  different 
fuel  conversion  rates  and  thus  different  durations  of  the  first  and  second 
phases  of  biomass  consumption.  The  variability  observed  may  be  relat¬ 
ed  to  operational  problems  arising  from  the  use  of  too  low  or  high  loads. 
Low  fuel  loads  have  contributed  to  lower  combustion  temperatures  and 
reduced  char  beds,  causing  difficulties  on  the  ignition  of  the  next  fuel 
batch  (combustion  chamber  recharge).  Thus,  it  is  generally  observed 
that  during  the  first  refuelling  after  heating  the  chamber  with  a  medium 
load  of  fuel,  temperatures  were  moderate  in  the  chamber,  ignition  oc¬ 
curred  easily,  and  the  cycle  presented  a  short  duration  (about  30 
min).  In  subsequent  batches,  ignition  was  delayed,  since  the  bed  of 
hot  coal  was  reduced,  and  the  cycle  took  twice  as  much  time.  However, 
after  performing  multiple  batches,  the  temperature  rose  whether  in  the 
combustion  chamber  or  in  the  char  bed.  For  high  fuel  loads,  the  combus¬ 
tion  cycles  lasted  from  40  to  80  min.  When  the  firebox  is  highly  loaded 
with  finely  split  pieces  of  wood,  excessive  compaction  may  occur.  Be¬ 
cause  the  combustion  air  cannot  reach  all  the  pieces,  full  ignition  may 
be  delayed.  On  the  other  hand,  in  the  tests  where  no  cleavage  was 
applied  to  the  logs,  difficulties  in  lighting  the  fire  were  again  felt  and 
low  fuel  consumption  rates  were  observed.  Larger  pieces  make  it  diffi¬ 
cult  to  ignite  because  they  have  larger  surface  areas,  making  them 
tougher  to  catch  fire. 

The  use  of  different  ignition  techniques  produced  no  significant 
effects  whether  on  the  combustion  air  inlet  or  on  the  biomass  consump¬ 
tion.  Since  the  load  used  to  start-up  the  fires  is  equivalent  to  those  used 
for  evaluating  the  emissions  of  medium  fuel  loads,  the  air  inlet  and  bio¬ 
mass  consumption  profiles  are  quite  similar  in  both  cases.  However,  it  is 
necessary  to  take  into  account  that  the  cold  start  involves  lower  temper¬ 
atures  and,  consequently,  lower  combustion  air  flows  at  the  beginning 
of  the  combustion  cycle.  Comparing  the  average  combustion  air  flow 
for  the  hot  start-up  (medium  load)  and  cold  start-up  experiments, 
only  slight  differences  were  observed  (about  1  Nm3  h”1  higher  when 
initiating  hot  start-up  combustion  tests). 

3.2.  Caseous  emissions 

The  emission  factors  were  calculated  following  the  methodology 
described  by  Calvo  et  al.  [28],  taking  into  account  the  concentration  of 
the  compound  in  the  flue  gas,  the  mean  combustion  gas  flow  rate  in 
the  chimney,  the  sampling  time  and  the  fuel  mass  consumed  (dry 


basis).  Fig.  5  displays  CO  and  THC  emission  factors,  for  pine  and  beech 
combustion,  obtained  for  each  condition  tested.  The  highest  CO  emis¬ 
sion  factor  (73.9  g  kg-1  of  wood  burned,  dry  basis)  for  pine  occurred 
when  the  combustion  at  high  load  conditions  with  non-split  logs  was 
carried  out.  Paired  t- tests  indicate  no  statistical  difference  (a  =  0.05) 
between  medium  load  and  either  high  (S  and  NS)  or  low  loads.  Howev¬ 
er,  the  difference  between  split  and  non-split  high  loads  was  statistically 
significant  ( p  =  0.037).  For  beech,  the  condition  that  produced  the 
highest  CO  emission  factor  (72.6  g  kg-1  of  wood  burned,  dry  basis) 
was  the  low  fuel  load.  The  difference  between  the  medium  load 
and  the  low  load  was  statistically  significant  (p  =  0.0233).  The  paired 
t- tests  performed  with  the  mean  CO  emission  factors  of  high  loads  (S 
and  NS)  indicated  no  statistical  difference  (a  =  0.05).  As  regards  the 
ignition  technique  (from  the  bottom  and  from  the  top),  for  beech  com¬ 
bustion,  the  p  value  revealed  that  there  is  a  significant  difference  be¬ 
tween  the  two  means  (p  =  0.0036).  Similarly  to  CO,  the  highest  THC 
emission  factor  for  pine  combustion  occurred  for  high  load  conditions 
without  cleavage  (9.0  g  kg~ 1  of  wood  burned,  dry  gases).  The  difference 
observed  between  medium  and  high  load  (NS)  is  statistically  significant 
(p  =  0.0063).  The  difference  observed  between  operation  with  medium 
loads  and  operation  with  secondary  air  supply  was  also  found  to  be 
statistically  significant  (p  =  0.0051).  For  beech  combustion,  the  condi¬ 
tion  that  generated  the  highest  THC  emission  factor  was  low  load 
(10.7  g  kg-1  of  wood  burned,  dry  gases).  The  difference  in  the  THC 
emission  factor  between  the  operation  with  medium  and  low  loads 
was  found  to  be  statistically  significant  (p  =  0.0012).  The  paired  t-tests 
performed  with  high  loads  (S  and  NS)  indicated  no  statistical  difference 
(a  =  0.05).  The  use  of  high  loads,  with  and  without  cleavage  of  the  logs, 
may  lead  to  incomplete  combustion.  This  occurs  because  a  large  fuel 
batch  increases  the  rate  of  devolatilisation  of  the  fuel,  which  causes  an 
insufficiency  on  the  combustion  air  supply  [8],  In  equipment  with  no 
secondary  combustion  chamber,  the  load  is  a  critical  aspect,  since  the 
empty  space  in  the  upper  part  of  the  chamber  acts  as  a  secondary  com¬ 
bustion  chamber.  When  using  an  excessively  high  fuel  load  this  space  is 
reduced  and  consequently  the  residence  time  is  not  enough  for  mixture 
of  the  flue  gas  with  combustion  air  and  gas  phase  oxidation  to  occur 
[42].  The  CO  and  THC  emission  factors  obtained  in  the  present  work 
for  beech  combustion  under  medium  load  conditions  are  in  general 
agreement  with  the  values  presented  by  Evtyugina  et  al.  [43]  for  the 
same  wood  fuel  and  combustion  appliance.  Amorim  et  al.  [44]  assessed, 
among  other  variables,  the  influence  of  biomass  size  and  combustion  air 
flow  rate  on  gaseous  emission  factors.  The  statistical  analysis  performed 
by  the  authors  revealed  that  the  most  influential  parameter  was  the 
biomass  size.  The  C02  emission  factors  were  higher  when  using  biomass 
with  smaller  diameter.  Conversely,  the  CO  and  CH4  emission  factors 
were  higher  when  combusting  with  larger  log  sizes.  This  trend  was 
also  observed  in  this  study.  The  emission  factors  presented  by  Amorim 
et  al.  [44]  are  close  to  the  values  obtained  in  the  present  study  for  beech 
combustion.  The  lower  emissions  of  incomplete  combustion  products 
when  using  smaller  pieces  result  from  the  fact  that  they  produce  a  lon¬ 
ger  flame  combustion  phase  and  a  shorter  smouldering  phase.  Hytonen 
&Jokiniemi  [40]  conducted  tests  in  a  combustion  stove  with  primary 
and  secondary  air  inlet.  The  tests  evaluated  the  influence  of  moisture 
content  (10  to  25%),  fuel  load  (1,  2  and  3  kg),  and  degree  of  cleavage. 
They  found  that  the  emissions  cannot  be  reduced,  even  by  modifying 
the  air  combustion  flow,  if  a  low  fuel  load  is  used.  In  these  conditions, 
the  combustion  intensity  was  low,  not  allowing  the  mixing  of  the  flue 
gases  and  combustion  air.  At  high  loads,  they  observed  high  combustion 
rates,  making  it  possible  to  obtain  an  adequate  mixing  degree,  although 
combustion  air  needs  have  been  high.  The  highest  degree  of  cleavage 
(over  30%  of  surface  area)  and  lower  moisture  content  generated  great¬ 
er  CO  emissions  for  high  fuel  loads.  When  applying  a  high  degree  of 
cleavage  at  low  loads  and  increased  moisture  content,  a  reduction  in 
emissions  was  observed.  The  explanation  lies  in  the  fact  that,  when 
the  reaction  area  is  high  and  the  moisture  content  is  much  reduced, 
the  combustion  is  more  intense  and  can  lead  to  local  conditions  of 
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Fig.  5.  Average  and  standard  deviations  of  CO  and  THC  emission  factors  (g  kg  1  wood  burned,  dry  basis). 


oxygen  deficit.  The  high  moisture  content  controls  the  combustion 
rates,  producing  lower  emissions  of  unburned  products. 

In  the  kind  of  appliance  used  in  the  present  work,  all  the  combustion 
air  is  supplied  as  primary  air  under  the  grate  of  the  fixed  bed.  The  sec¬ 
ondary  inlet  allows  the  supply  of  secondary  air  over  the  firewood 
batch  and,  consequently,  an  improved  combustion  of  the  devolatilised 
gases,  leading  to  lower  emissions  compared  to  those  of  the  appliance 
with  the  traditional  design. 

Only  a  slight  reduction  in  the  THC  emission  factor  was  observed 
when  using  the  secondary  air  inlet.  This  may  have  to  do  with  the  fact 
of  not  having  done  preheating  of  the  combustion  air  and  not  having 
used  secondary  air  jets.  The  secondary  combustion  air  preheating  pre¬ 
vents  the  temperature  reduction  in  the  combustion  chamber  and  the 
secondary  air  jets  enhance  the  mixing  of  devolatilised  gases  and  com¬ 
bustion  air  [41],  However,  it  is  necessary  to  take  into  account  that 
these  tests,  rather  than  to  design  an  optimised  stove,  intended  to  simu¬ 
late  a  low  cost,  easy  to  implement,  measure  to  reduce  emissions  from 
residential  wood  combustion. 

Some  variability  in  emission  factors  was  observed,  which  may  derive 
from  several  factors.  There  is  variation  in  the  biomass  composition,  even 
within  the  same  species,  which  can  induce  differences  in  gaseous  emis¬ 
sion  factors.  Combustion  conditions  and  operation  are  also  very  impor¬ 
tant  factors.  For  example,  high  rates  of  combustion  result  in  increased 
emissions  of  unburned  products,  which  can  be  induced  by  the  low 
moisture  content  of  the  biofuel  [17].  The  distribution  of  biomass  itself 
on  the  grate  and  degree  of  compaction  can  cause  significant  differences 
in  emissions  [45  j. 

The  C0/C02  ratio  is  a  measure  of  the  relative  combustion  efficiency 
in  terms  of  conversion  of  the  fuel.  A  high  ratio  reflects  lower  combustion 
efficiency.  The  typical  ratio  for  flame  combustion  phase  is  less  than  0.1 
[29].  On  average,  the  ratios  obtained  for  all  combustion  experiments 
were  below  this  value,  which  indicates  that  the  flame  combustion 
phase  was  dominant.  The  highest  ratios  were  obtained  for  the  tests 
with  low  load  operation  for  both  biofuels. 

There  is  a  clear  relationship  between  CO  and  C02  emission  factors 
and  combustion  temperatures.  Thus,  CO  emission  factors  increase 
with  decreasing  combustion  temperatures,  whilst  the  C02  emission 


factors  have  an  opposite  behaviour,  increasing  with  increasing  combus¬ 
tion  temperatures.  The  technology  of  combustion  is  therefore  a  key  fac¬ 
tor  in  the  emissions  of  these  compounds.  Advanced  combustion 
technologies  allow  the  reduction  of  CO,  enhancing  its  oxidation  to 
C02.  On  the  other  hand,  fireplaces,  where  there  is  a  high  combustion 
air  inlet  reducing  combustion  temperatures,  originate  high  CO  emission 
factors. 

3.3.  PM10  emission  factors 

The  emission  factors  were  calculated  following  the  methodology  de¬ 
scribed  by  Calvo  et  al.  [28],  taking  into  account  the  concentration  of 
PM10  in  the  dilution  tunnel,  the  mean  combustion  gas  flow  rate  also  in 
the  dilution  tunnel,  the  sampling  time  and  the  fuel  mass  consumed 
(dry  basis).  The  emission  factors  for  “softwood"  combustion  ranged  be¬ 
tween  5.07  ±  0.30  and  16.74  ±  1.92  (gkg-1  of  wood  burned,  dry  basis). 
With  regard  to  the  load  and  degree  of  cleavage,  the  lowest  PM10  emis¬ 
sion  factor  was  observed  when  operating  the  stove  with  medium  load, 
whilst  the  highest  was  generated  when  loading  a  lower  amount  of 
wood.  The  difference  found  was  more  than  twice,  as  can  be  seen  in 
Fig.  6.  This  difference  is  statistically  significant  (p  <  0.0001).  The  t-test 
revealed  that  the  difference  found  between  the  mean  PM10  emission 
factors  for  medium  and  high  loads  (S  and  NS)  was  statistically  signifi¬ 
cant  (p  =  0.0121  and  p  =  0.0004,  respectively).  Too  small  chips  and 
logs  should  only  be  used  for  the  ignition.  When  small  logs  are  used  to 
feed  the  stove,  the  high  devolatilisation  rate  leads  to  oxygen  deficiency, 
causing  high  pollutant  emissions  [46],  Overall,  non-split  wood  gener¬ 
ates  higher  emissions  than  split  wood.  The  start-up  emissions  can  be 
substantial,  since  the  temperatures  in  combustion  chamber  are  reduced 
and  the  combustion  rate  is  low,  leading  to  a  high  PM  emission  factor. 
The  duration  of  this  step  depends  largely  on  the  fuel  properties,  such 
as  moisture  content,  density  and  size  [47],  It  could  be  verified  that  the 
use  of  big  logs  extends  the  initial  combustion  stage,  contributing  to 
higher  PM  emissions.  The  top-down  ignition  can  reduce  the  PM10  emis¬ 
sion  factor  to  less  than  half  when  compared  with  the  traditional  tech¬ 
nique.  The  differences  between  the  mean  PM10  emission  factors  for 
both  lighting  techniques  are  statistically  significant  (p  =  0.0016). 
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Fig.  6.  Average  and  standard  deviations  of  PM10  emission  factors. 


However,  the  difference  between  the  mean  PM10  emission  factors  for 
medium  load  (hot  start)  and  top  ignition  (cold  start)  is  not  statistically 
significant  (p  =  0.7986).  On  the  other  hand,  when  the  mean  PM]0  emis¬ 
sion  factor  for  medium  load  is  compared  with  that  of  the  bottom  igni¬ 
tion  technique  the  difference  has  shown  to  be  statistically  significant 
(p  <  0.0001 ).  The  traditional  ignition  technique  from  the  bottom  causes 
the  simultaneous  firing  of  the  entire  batch  of  fuel,  leading  to  high  com¬ 
bustion  rates  and  therefore  oxygen  deficit  zones  in  the  combustion 
chamber,  which  results  in  incomplete  combustion.  On  the  other  hand, 
the  top-down  ignition  promotes  a  gradual  combustion  of  the  fuel 
batch  resulting  in  more  complete  combustion  [48].  The  use  of  a  second¬ 
ary  combustion  air  flow  was  the  operating  condition  that  allowed 
obtaining  the  lowest  PM10  emission  factor.  The  soot  and  condensable 
organic  matter  can  be  reduced  by  staged  combustion  using  primary 
air  and  a  subsequent  gas  phase  oxidation  with  secondary  air  [49]. 

PM]0  emission  factors  obtained  for  beech  combustion  for  each  oper¬ 
ating  condition  are  presented  in  Fig.  6.  The  emission  factors  for  this 
hardwood  ranged  from  4.86  ±  1.06  to  11.0  ±  2.50  (g  kg"1  of  wood 
burned,  dry  basis)  for  the  different  loads  tested.  Contrary  to  what  was 
observed  for  the  softwood  combustion,  it  seems  that,  except  for  low 
loads,  the  other  operating  conditions  present  similar  emission  factors. 
The  differences  between  the  means  of  PM10  emission  factor  for  medium 
load  were  statistically  significant  when  compared  with  the  operation 
with  low  (p  =  0.0300)  and  high  load  (S)  (p  =  0.0033).  Compared  to 
bottom  lighting  (6.06  ±  1.65  g  kg"1  of  wood  burned,  dry  basis),  a  slight 
reduction  of  the  emission  factor  was  observed  when  the  ignition  of  the 
hardwood  was  carried  out  from  the  top  (5.24  ±  1.10  g  kg"1  of  wood 
burned,  dry  basis),  but  this  difference  was  not  statistically  significant 
(p  =  0.4938).  As  observed  for  the  combustion  of  pine,  the  highest 
emission  factors  were  generated  under  low  load  conditions.  However, 
contrary  to  what  was  registered  for  the  softwood  combustion,  the  low¬ 
est  emission  factor  was  generated  when  operating  under  high  loads 
without  cleavage.  The  lower  emission  factor  may  be  associated  with  a 
higher  density  of  hardwood  that  promotes  lower  combustion  rates 
compared  to  those  observed  for  softwood.  Hardwood  burns  slower 
and,  thus,  prevents  the  occurrence  of  oxygen  deficit  zones  in  the  com¬ 
bustion  chamber.  The  high  temperatures  recorded  during  operation 


with  high  loads,  associated  with  good  combustion  air  supply,  promote 
a  more  efficient  combustion  [45].  Similarly  to  what  was  observed  in 
this  study  for  the  hardwood  tests,  Schmidl  et  al.  [24]  found  that  combus¬ 
tion  with  “high  fuel  load"  in  a  sophisticated  stove  presented  lower  PM10 
emissions  compared  with  the  “normal”  operation  condition.  The  au¬ 
thors  point  out  the  reduced  air  combustion  flow  as  a  critical  factor  in 
particulate  matter  emissions.  Some  authors  have  already  tried  to  assess 
the  influence  of  the  fuel  load  and  oxygen  supply  on  particulate  emis¬ 
sions.  Ozgen  et  al.  [50]  tested  three  closed  manual  appliances,  namely 
a  closed  fireplace,  a  traditional  stove,  and  an  advanced  stove.  These  ap¬ 
pliances  were  tested  to  evaluate  the  influence  of  increased  fuel  load  and 
decreased  air  through  the  application  of  two  cycles  (A  and  B).  Both  cy¬ 
cles  comprised  an  ignition  phase  and  three  consecutive  fuel  loadings. 
The  difference  between  cycles  derived  from  the  fact  that,  together 
with  the  complete  closure  of  the  primary  air  valve,  the  last  fuel 
reloading  in  A  was  1.5  times  higher  than  the  nominal  value  in  B.  For 
the  advanced  stove,  the  authors  observed  a  marked  difference  in  the 
PM  emissions  between  the  two  cycles.  The  difference  in  the  average 
PM  emissions  was  not  so  noticeable  for  the  tests  with  the  traditional 
stove  and  the  closed  fireplace.  However,  median  values  lower  than 
mean  values  indicated  the  presence  of  very  high  emission  factors  during 
some  test  runs.  In  accordance  with  these  results,  Shen  et  al.  [51  ]  ob¬ 
served  that  under  different  feeding  rates,  PM,  OC  and  EC  emission  fac¬ 
tors  were  not  statistically  different.  However,  they  noticed  that  PAH 
emission  factors  and  their  oxygenated  and  nitrated  derivatives  were 
much  higher  in  the  fast  burning. 

Tissari  et  al.  [21]  conducted  tests  on  a  conventional  masonry  heater 
in  which  two  conditions  have  been  tested,  namely  “normal  combus¬ 
tion”  and  “smouldering  combustion”.  The  “smouldering  combustion” 
conditions  have  been  created  using  a  high  fuel  load  and  by  restricting 
the  entry  of  combustion  air.  For  normal  combustion  conditions  the  au¬ 
thors  obtained  an  emission  factor  of  1.8  ±  0.5  PMi  g  kg"1  of  wood 
burned.  For  “smouldering  combustion”  conditions,  the  emission  factor 
was  about  six  times  higher,  with  a  mean  value  of  11.1  ±  3.9  g  kg"1. 
This  value  is  identical  to  that  obtained  in  this  study  under  low  load  con¬ 
dition  for  beech  combustion  (1 1.0  g  kg"1  of  wood  burned,  dry  basis). 
Thus,  reduced  combustion  temperatures,  poor  mixture  between  the 
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flue  gas  and  combustion  air  or  short  residence  times  contributed  to  in¬ 
creased  emissions  of  particulate  matter.  The  particle  emission  factors 
found  in  the  literature  may  show  great  variability,  due  to  the  fuel  char¬ 
acteristics,  combustion  equipment  design,  combustion  conditions,  dilu¬ 
tion  techniques  and  sampling  procedures  [21,22,32,42,52], 

3.4.  Number  emission  factors 

The  global  number  emission  factors  for  each  operating  condition  and 
for  both  fuels  are  presented  in  Table  3.  With  regard  to  the  combustion 
of  pine  with  different  batch  sizes,  the  maximum  global  emission  factor 
was  recorded  when  operating  with  low  load  (1.05  x  1017  particles -kg^1 
biomass),  whereas  the  minimum  was  registered  when  burning  non-split 
logs  in  high  load  conditions  (4.08  x  1016  ±  3.48  x  1015  particles -kg^1  of 
wood  burned).  The  operation  with  secondary  air  originated  a  low 
emission  factor  (2.19  x  1016  ±  1.50  x  1015  particles- kg-1  of  wood 
burned)  compared  with  the  emission  factors  mentioned  above. 
Beech  combustion  showed  the  same  trend  described  above.  The 
overall  maximum  emission  factor  was  observed  when  operating 
with  low  load  (1.09  x  1017  ±  1.08  x  1016  particles- kg-1  of  wood 
burned),  whereas  the  minimum  occurred  under  conditions  of  high 
load  (S)  (2.47  x  1016±  6.09  x  1014  particles- kg-1  ofwood  burned). 

Regarding  the  lightning  technique,  it  can  be  seen  that,  for  both  woods, 
the  overall  emission  factor  is  higher  when  the  ignition  is  carried  out  from 
the  top,  although  in  the  case  of  beech  the  difference  between  techniques 
is  negligible.  The  top-down  ignition  technique  generated  emission  factors 
between  5.26  x  1016  ±  1.15  x  1016  and  3.68  x  1016  ±  2.54  x  1015 
(particles -kg-1  of  wood  burned)  for  pine  and  beech,  respectively. 
The  bottom  ignition  technique  contributed  to  emission  factors  between 
5.72  x  1015  and  3.69  x  1016  ±  2.33  x  1015  (particles -kg”1  of  wood 
burned)  for  pine  and  beech,  respectively.  The  top  ignition  leads  to  a  pro¬ 
gressive  combustion  of  the  fuel,  resulting  in  more  complete  combustion, 
as  described  above.  Thus,  the  increased  combustion  efficiency  leads  to 
number  emission  factors  exceeding  the  number  recorded  for  the  tradi¬ 
tional  lighting  technique,  which  originates  incomplete  combustion  condi¬ 
tions.  As  observed  for  the  PM  emission  factors,  this  trend  is  remarkably 
visible  for  the  combustion  of  pine.  The  increase  in  the  combustion  tem¬ 
perature  leads  to  the  emission  of  a  larger  number  of  particles. 

Tissari  [8]  observed  that  when  the  batch  size  was  doubled  in  a  con¬ 
ventional  masonry  heater,  all  emissions  increased  except  the  particle 
number  emission.  In  our  study,  the  same  trend  was  verified  in  the  com¬ 
bustion  of  pine.  Using  identical  batches  and  modifying  the  size  of  the 
logs,  Tissari  [8]  observed  that  the  particle  number  emissions  from 
small  logs  were  1.4-fold  higher  than  from  big  logs.  In  the  present 
study,  for  the  same  batch  sizes,  the  number  emission  factor  was  higher 
when  using  logs  without  cleavage  compared  to  the  emission  factor  ob¬ 
tained  with  split  logs. 

Tissari  et  al.  [22]  studied  the  influence  of  the  type  of  equipment  and 
operating  conditions  on  the  particle  number  emission  factors.  The  exper¬ 
iments  involved  the  use  of  combustion  appliances  typically  employed  in 
Finland  (conventional  masonry  stove  and  modern  stove).  Different  bio¬ 
mass  loads  and  log  sizes  were  evaluated.  The  authors  concluded  that 
the  number  emission  factor  was  higher  for  the  modern  combustion 


equipment  (5.9  x  1014  particles -kg”1  ofwood  burned)  than  that  obtain¬ 
ed  with  the  conventional  stove  (3.1  x  1014  particles  •  kg“ 1  of  wood 
burned).  Moreover,  the  authors  found  that  higher  loads  produce  smaller 
number  of  particles  (2.0  x  10 14  particles  -kg-1  biomass)  in  comparison 
with  low  loads  (4.0  x  1014  particles -kg-1  of  wood  burned).  This  trend 
was  also  observed  in  the  present  study,  although  the  recorded  values 
are  higher.  The  variations  are  probably  related  to  differences  in  character¬ 
istics  of  combustion  devices.  When  the  degree  of  cleavage  applied  was 
higher,  i.e.  when  smaller  wood  pieces  were  used,  the  authors  found  an  in¬ 
crease  in  the  number  of  particles  (3.1  x  1014  particles -kg-1  of  wood 
burned)  comparatively  to  the  number  of  particles  emitted  during  the 
combustion  of  larger  pieces  of  wood  (2.2  x  1014  particles -kg^1  of  wood 
burned)  [20]. 

Tissari  et  al.  [21]  conducted  tests  on  a  conventional  masonry  stove. 
These  tests  comprised  “normal  combustion”  and  “smouldering  combus¬ 
tion"  conditions.  The  slow  combustion  conditions  have  been  created 
using  an  excess  of  fuel  load  and  restricting  the  entry  of  combustion 
air.  The  authors  found  that  in  improper  combustion  conditions  (i.e. 
with  excess  of  fuel  and  oxygen  deficit  in  the  combustion  chamber), 
the  number  of  particles  was  about  20%  lower  compared  to  "normal” 
combustion  conditions.  The  smaller  number  of  particles  occurs  due  to 
coagulation  and  condensation  that  contribute  to  the  growth  of  the  parti¬ 
cles.  Thus,  the  particle  number  is  reduced  and,  moreover,  size  increases. 
The  authors  found  that  the  mass  emission  was  quite  high  due  to  the 
large  sizes  of  particles  emitted.  Wardoyo  [53]  investigated  the  emissions 
of  particles  in  a  modified  stove  in  order  to  allow  the  introduction  of  a  con¬ 
trolled  amount  of  air  and  variable  combustion  rates.  The  author  tested 
“vigorous  combustion"  and  “smouldering  combustion”  conditions.  For 
the  first  condition  a  linear  relationship  (R2  >  0.93)  between  temperature 
and  particle  number  concentration  was  found.  The  particle  number  for 
the  “vigorous  combustion"  conditions  ranged  from  3.3  to  5.7  x  1015 
( particle  •  kg~ 1  of  wood  burned ).  The  “smouldering”  conditions  generated 
a  lower  number  of  particles  (from  2.8  to  44.8  x  1013  particles -kg-1  of 
wood  burned),  as  observed  by  Tissari  et  al.  [21  ].  The  tests  performed  in 
this  study  produced  fewer  particles  for  high  load  conditions  with  non¬ 
split  logs,  for  which  the  combustion  air  supply  is  insufficient  for  high 
devolatilisation  rates  of  biomass.  Note  that,  for  operation  under  low 
load,  the  PM  mass  emission  factors  were  the  highest  among  the  tested 
conditions.  As  expected,  these  results  were  associated  with  unfavourable 
combustion  conditions  with  reduced  temperatures  and  low  combustion 
rates.  It  should  be  taken  into  account,  however,  that  the  exact  replication 
of  conditions  between  successive  combustion  experiments  is  difficult, 
which  causes  variations  in  results. 


3.5.  Aerosol  size  distributions 

Regarding  the  fuel  load  during  pine  combustion,  the  highest  dg  was 
registered  in  the  operation  with  high  load  (NS)  (0.265  ±  0.008  pm).  The 
operation  with  secondary  air  supply  generated  the  lowest  Dg  (0.153  ± 
0.000  pm)  for  the  combustion  of  the  coniferous  wood.  For  beech  com¬ 
bustion,  the  highest  Dg  was  observed  when  operating  the  stove  with 
high  fuel  load  (S)  (0.243  ±  0.046  pm). 


Table  3 

Mean  geometric  diameter  (Dg)  and  number  emission  factors  ( particles  ■  kg  1  ofwood  burned)  for  each  operating  condition  (mean  ±  standard  deviation)  for  both  fuels. 


Fagus  sylvatica  Pinus  pinaster 


Dg(pm)  EF  (particles -kg  1  ofwood  burned)  Dg  (pm)  EF  (particles -kg  1  of  wood  burned) 


Low  load 

0.205 

± 

0.016 

1.09E 

+ 

17 

± 

1.08E 

+ 

16 

0.175 

± 

NA 

1.05E 

+ 

17 

± 

NA 

Medium  load 

0.234 

± 

0.011 

3.75E 

+ 

16 

± 

9.87E 

+ 

14 

0.175 

± 

0.037 

4.08E 

+ 

16 

± 

1.29E 

+ 

15 

High  load  (S) 

0.243 

± 

0.046 

2.47E 

+ 

16 

± 

6.09E 

+ 

14 

NA 

NA 

High  load  (NS) 

0.227 

± 

0.020 

3.09E 

+ 

16 

± 

5.96E 

+ 

15 

0.265 

± 

0.008 

4.08E 

+ 

16 

± 

3.48E 

+ 

15 

Secondary  combustion  air  supply 

NA 

NA 

0.153 

± 

0,000 

2.19E 

+ 

16 

± 

1.50E 

+ 

15 

Top  ignition 

0.165 

± 

0.001 

3.68E 

+ 

16 

± 

2.54E 

+ 

15 

0.149 

± 

0.016 

5.26E 

+ 

16 

± 

1.15E 

+ 

16 

Bottom  ignition 

0.178 

± 

0.022 

3.69E 

+ 

16 

± 

2.33E 

+ 

15 

0.158 

± 

NA 

5.72E 

+ 

15 

± 

NA 

NA  -  not  applicable. 
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With  regard  to  the  ignition  technique,  it  has  been  found  that,  for 
pine  combustion,  the  lowest  Dg  was  recorded  when  using  the  top  igni¬ 
tion  technique  (0.149  ±  0.016  pm).  For  beech  combustion,  the  lowest 
Dg  was  observed  when  the  ignition  was  performed  from  the  top 
(0.165  ±  0.001  pm). 

Particle  size  increases  during  conditions  of  incomplete  combus¬ 
tion  due  to  condensation  and  agglomeration  processes  induced  by 
the  low  temperatures  and  poor  mixing.  Tissari  et  al.  [22]  found  that 
modern  biomass  combustion  equipment  emits  more  particles,  but 
smaller  diameters.  They  also  observed  that  the  use  of  high  biomass 
loads  causes  the  increase  of  the  particle  diameter.  This  trend  was 
also  observed  when  using  small  pieces  of  wood.  Pettersson  et  al. 
[17]  concluded  that  the  particulate  material  resulting  from  domestic 
combustion  processes  is  dominated  by  submicron  particles  in  the 
range  from  0.1  to  0.3  pm.  The  researchers  also  found  differences  in 
the  size  distributions  for  the  distinct  stages  of  the  combustion 
cycle.  In  this  work,  these  trends  were  also  observed.  Conditions  of  in¬ 
complete  combustion  cause  the  formation  of  larger  aggregates  of 
soot  and  condensation  of  organic  material  into  the  pre-existing 
particles  during  the  cooling  of  the  exhaust  effluent,  modifying  the 
particle  diameters  [17].  The  increase  of  particle  diameters  under 
conditions  of  incomplete  combustion  was  also  observed  by  other  au¬ 
thors.  Tissari  et  al.  [21]  described  a  2.5-fold  increase  in  the  geometric 
mean  diameter  for  incomplete  combustion  conditions  when  com¬ 
paring  with  "normal”  combustion  conditions. 

The  study  of  the  temporal  evolution  of  the  number  of  particles  in 
each  channel  has  also  been  done.  Analysing  the  operation  with  medium 
load  (Fig.  7)  it  was  observed  a  maximum  number  of  9.6  x  106 
particles -cm-3  for  the  size  bin  of  0.28  pm  during  the  first  stage  of  the 
combustion  cycle.  Subsequently,  there  was  a  progressive  decrease  in 
the  number  of  particles  collected  on  the  following  channels.  During 
the  flaming  stage,  when  the  number  of  particles  is  at  its  maximum,  a 
peak  of  1.9  x  107  particles-cm~3  for  the  size  bin  of  0.22  pm  was  regis¬ 
tered.  During  the  flameless  combustion  stage,  the  peak,  7.1  x  107 

a) 


Diameter  (pm) 


b) 


Diameter  (pm) 

Fig.  7.  Typical  evolution  of  the  distribution  of  the  aerosol  from  the  heating  up  combustion 
stage  (stage  1 )  until  the  flameless  combustion  stage  (stage  3)  during  the  combustion  of 
(a)  beech  combustion  operating  with  medium  load  (b)  pine  combustion  operating  with 
medium  load. 


particles -cm-3,  was  recorded  for  the  first  channel  (0.08  pm).  At  this 
stage,  the  number  of  particles  within  the  smallest  size  ranges  experi¬ 
enced  an  increase  compared  with  the  values  observed  for  phases  1 
and  2.  In  the  case  of  beech,  the  evolution  of  the  number  of  particles  is 
similar  to  that  described  above  for  pine  combustion.  In  the  first  stage 
of  the  combustion  cycle,  a  peak  of  8.1  x  106  particles- cm-3  in  the  size 
channel  of  0.31  pm  was  recorded.  Thereafter,  there  is  a  progressive  re¬ 
duction  in  the  number  of  particles.  During  the  flaming  phase,  a  maxi¬ 
mum  of  1.3  x  107  particles -cm- 3  in  the  size  channel  of  0.25  pm  was 
achieved.  For  flameless  combustion,  the  largest  number,  6.5  x  106 
particles -cm-3,  was  recorded  for  the  third  channel  (0.13  pm).  At  this 
stage,  an  increase  in  the  number  of  particles  within  the  smallest  size 
bins  is  observed  compared  to  the  values  of  phases  1  and  2. 


4.  Conclusions 

The  present  work  had  as  objective  to  experimentally  characterise 
and  quantify  the  gaseous  and  particulate  emissions  during  combustion 
in  a  residential  woodstove  with  variations  in  fuel  and  operating  condi¬ 
tions.  The  comparison  of  values  from  this  study  with  literature  data 
showed  differences,  depending  mainly  on  operating  procedures  and 
combustion  technologies.  Thus,  to  implement  measures  and  as  a  guide¬ 
line  for  future  regulations,  emission  factors  from  different  types  of  com¬ 
bustion  devices  under  different  operation  modes  are  needed.  This  work 
may  be  potentially  useful  because  it  provides  some  guidance  on  best 
practices  for  households  and  end-users  in  Southern  European  countries, 
where  the  combustion  of  wood  in  traditional  appliances  is  becoming 
ever  more  common. 

The  analysis  of  the  CO  and  THC  evolution  over  a  combustion  cycle 
showed  that  the  initial  stage  is  critical  from  the  point  of  view  of  un¬ 
burnt  compound  emissions,  due  to  the  rapid  devolatilisation  of  bio¬ 
mass  in  a  short  period  of  time.  The  lowest  particle  emissions  were 
generated  under  operation  with  medium  fuel  loads,  whereas  the 
highest  were  achieved  for  low  and  high  loads  of  pine  and  beech,  re¬ 
spectively.  The  secondary  air  supply  enables  to  slightly  reduce  the 
THC  emission  factor.  Cleavage  also  demonstrated  to  play  an  impor¬ 
tant  role.  Particulate  emissions  increase  with  the  decrease  of  the  sur¬ 
face  area  of  the  logs.  Regarding  the  ignition  technique,  it  was 
observed  that  the  top-down  ignition  technique  can  reduce  the 
PM10  emission  factor  to  less  than  half  compared  with  the  traditional 
technique.  The  highest  particle  number  emission  factors  were  gener¬ 
ated  for  the  combustion  of  low  load  of  both  fuels.  The  lowest  number 
emission  factors  were  observed  when  operating  with  high  load  (S), 
for  the  combustion  of  beech,  and  when  using  a  secondary  air  supply, 
for  the  combustion  of  pine.  The  smallest  global  particle  geometric 
mean  diameters  were  reached  under  low  fuel  load  operation  for 
beech  and  during  the  use  of  secondary  air  supply  for  pine  combus¬ 
tion.  The  use  of  high  loads,  with  split  and  non-split  logs,  generated, 
for  both  woods,  the  highest  global  Dg.  Regarding  the  ignition  tech¬ 
nique,  the  global  geometric  mean  diameters  were  the  lowest 
among  all  the  conditions  tested  using  both  techniques. 

In  batch  operation,  the  conditions  are  highly  variable  throughout  the 
different  combustion  phases.  The  log  size  affected  emissions  more  than 
did  batch  size.  The  present  study  showed  that  a  large  variation  in  emis¬ 
sions  can  easily  be  achieved  by  the  woodstove  operators.  Among  all  the 
conditions  tested,  it  was  possible  to  conclude  that  the  ignition  technique 
and  combustion  air  supply  are  the  operational  parameters  that  most  no¬ 
ticeably  affect  PM10  emissions. 
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